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O ■ Abstract 

°. 

CN \ In the framework of analytic approach to QCD the nonperturbative contributions 

£>Y in running coupling of strong interaction up to 4-loop order are obtained in an 

explicit form. For all Q > A they are shown to be represented in the form of 
an expansion in inverse powers of Euclidean momentum squared. The expansion 
coefficients are calculated for different numbers of active quark flavors n/ and for 
different number of loops taken into account. On basis of the stated expansion 
the effective method for precise calculation of the analytic running coupling can be 
' developed. 
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Q_i! It is widely believed that nonphysical singularities of the perturbation theory in the 
o ' infrared region of QCD should be canceled by the nonperturbative contributions. The 
nonperturbative contributions arise quite naturally in an analytic approach jl]] — || to 
QCD. The idea of the approach goes back to Refs. || [7[] devoted to the nonphysical ghost 
pole problem in QED. In recent papers it is suggested to solve the ghost pole problem 
in QCD demanding the running coupling constant be analytic in Q 2 (Q 2 is the Euclidean 
momentum squared). As a result of the procedure, instead of the one-loop expression 
a^(Q 2 ) = (47r/6 )/ ln(Q 2 /A 2 ) taking into account the leading logarithms and having the 
ghost pole at Q 2 = A 2 , one obtains new expression 



A 2 

+ 



ln(Q 2 /A 2 ) A 2 -Q 



Eq. ([!]) is an analytic function in the complex Q 2 -plane with a cut along the negative real 
semiaxis. The pole of the perturbative running coupling at Q 2 = A 2 is canceled by the 
nonperturbative contribution (A 2 ~ /i 2 exp{— 4ti / (b a an (fj, 2 ))} at a an (fi 2 ) — > 0) and the 
value a^(0) = 47r/6 appeared to be finite and independent of the normalization condi- 
tions (independent of A). The most important feature of the procedure is the stability 
property 0, [| of the value of the "analytically improved" running coupling constant at 
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zero with respect to high corrections, a^(0) = «^(0) = a^(0). This property provides 
the high corrections stability of a an (Q 2 ) in the whole infrared region. We shall demon- 
strate that a an (0) = 4n/b for any finite loop order initial a s (Q 2 ) in the standard inverse 
powers of logarithms expansion form being a consequence of the asymptotic freedom 
property. 

Note that for the one-loop case we have not only the spectral representation for a an (Q 2 ) 
but, in the first place, the nonperturbative contribution is extracted from the analytic 
running coupling explicitly, 

a an (Q 2 ) = a pt (Q 2 ) + a n *(Q 2 ), (2) 

and, in the second place, the 1-loop order nonperturbative contribution can be presented 
as convergent at Q 2 > A 2 of constant signs series in the inverse powers of the momentum 
squared, 

4tt^ {A 2 ^ n 



a 



where c n = —1. For the standard as well as for the iterative 2-loop perturbative input 
the nonperturbative contributions in the analytic running coupling are calculated explic- 
itly in Ref. ||. In the ultraviolet region the nonperturbative contributions also can be 
represented as a series in inverse powers of the momentum squared. The 3-loop case is 
considered in Ref. || where, to handle the singularities originating from the perturbative 
input, the method which is more general then that of Ref. |S| was developed. 

In this paper we extract in an explicit form the nonperturbative contributions to 
otaniQ 2 ) up to 4-loop order in analytic approach to QCD, and find the coefficients c n of 
the expansion of the form (|3[). This gives the effective method for a precise calculation of 
a an (Q 2 ) which is not connected with numerical integration. 

The behavior of the QCD running coupling a s (Q 2 ) is defined by the renormalization 
group equation 

Q2 ^p = P{as) = ^ + Pia * + + ^ + 0(c ^' (4) 

where the coefficients ]TI| — 



1 2 
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Here n/ is number of active quark flavors and ( is Riemann zeta- function, £3 = ((3) = 
1.202056903... . The first two coefficients /3q, (3% do not depend on the renormalization 
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scheme choice. The next coefficients do depend on this choice. Being calculated within 
the MS'-scheme in an arbitrary covariant gauge for the gluon field they appeared to be 
independent of the gauge parameter choice. Integrating Eq. (H) yields 



i lna,(Q 2 ) + ^ (fob ~ Pi) a s (Q 2 ) + ^ (ft 3 - 2/W2 + ffifc) « 2 (Q 



a s (Q 2 ) ^ ft— ' ' Pi 

+ O(a 3 s (Q 2 )) = -P ln(Q 2 /A 2 



+ C. 



(6) 



The integration constant is represented here as a combination of two constants A and 
C. Dimensional constant A is a parameter which defines the scale of Q and is used 



for developing the iteration procedure. Iteratively solving Eq. 
ln(Q 2 /A 2 ) — > 00 one obtains 



for a s {Q 2 ) at L 



PtL 2 



(InL + C) 2 - (InL + C) - 1 + 



P0P2 
PI 



+ 



Jl_ 

P$L* 



(InL + C) 3 - -(InL + C) 5 



W0P2 



^ PI )Q» L + W + 2 2PI 

( 7 ) 

where C = ln(— p ) + (P /Pi)C. Within the conventional definition of A as A^- [TJ]] 
one chooses C — 0. At that the functional form of the approximate solution for a s (Q 2 ) 
turns out to be somewhat simpler, but it requires distinct for different rif. With 
this choice Eq. (^) at three loop level corresponds to the standard solution written in 



/Soft 



+ " ' 



the form of the expansion in inverse powers of logarithms and at four loop level it 
corresponds to [16|. We shall deal with nonzero C because this freedom can be useful for 



an optimization of the finite order perturbation calculations. Moreover, in the presence of 
the n/-dependent constant C it is possible to construct matched solution of Eq. @ with 
universal rif independent constant A |Hfl . 

Let us introduce the function <&(z) of the form 



i_ h Hz) + c + p 



(ln(z) + C) 2 _ ln(z) + C k 



(Hz) + cf h {Hz) + cf . M£ + c_ 



+ 



where the coefficients b, k, and k are equal to 

Pi 



K 



PI 
-1 + 



h 2 ' 
°0 



26i 

A 2 



■1 + 



60&2 
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(9) 



To choose the main branch of the multivalued function (H) we cut complex z-plane along 
the negative semiaxis. Then the solution (|7|) can be written as a s (Q 2 ) = (4ft/bo)a(x), 
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where a(x) = <3>(lnx), x = Q 2 /A 2 . Function a(x) is unambiguously defined in the complex 
x-plain with two cuts along the real axis, physical cut from minus infinity to zero and 
nonphysical one from zero to unity. At x ~ 1 the perturbative running coupling has 
singularities of a different analytical structure. Namely, at x ~ 1 the leading singularities 
are 



a w (x) 



a®(x) 



x — 1 
b 2 



x-i: 



ln(x 



(x - iy 



ln 2 (x- V, 



b 3 



x-1) 



ln 3 (x - 1). 



(10) 



This is not an obstacle for the analytic approach which removes all this nonphysical 
singularities. 

The analytic running coupling is defined by the spectral representation 



da 



X 



a 



■Pi?) 



with the spectral density p(a) = Im a an (—cr — iO) = Ima(— a — iO) where a(x) is the per- 
turbative running coupling. It is seen that dispersively modified coupling of the form ( |TT]) 
has analytical structure which is consistent with causality. 

By making the analytic continuation of Eq. (|8|) into the Minkowski space by formal 
substitution x = —a — iO one obtains the spectral density as an imaginary part of $(lna — 
in). Function a(x) is regular and real for real x > 1. So, to find the spectral density p(c) 
we can use the reflection principle (a(x))* = a(x*) where x is considered as a complex 
variable. Then 



2i 



(<3>(ln<7 - in) - $(lncr + in)) . 



(12) 



By the change of variable of the form a = exp(t), the analytical expression is derived 
from (0), (pD as follows: 



Chan. \ X ) 



/ dt 

2ni 



x + e l 



x {$(t - in) - $(t + in)} . 



(13) 



Let us prove that a an (0) = 1. It follows from Eq. ( p~3| ) that 

1 °° 

a an (0) = — J dt {$(t - iir) - $(t + m)} 



1 

2ni 



dt 



m) — 



2ni 
1 

t — in 



$(t + in) 



1 



t + in 



+ 



1 



1 



t — in t + in 



(14) 



For the first term in Eq. (|I4[) close the integration contour in the lower half-plane of the 
complex variable t by the arch of the "infinite" radius without affecting the value of the 
integral. We can do it because the integrand multiplied by t goes to zero at | t |— > oo. 
There are no singularities inside the contour, so we obtain zero contribution from the term 
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considered. For the second term we close the integration contour in the upper half-plane 
of the complex variable t with the same result. So we have 



1 1 



dt 



t — iix t + m 



1. (15) 



For any finite loop order the expansion structure of the perturbative solution in inverse 
powers of logarithms ensure the property of the analytical coupling a an (0) = 1. The 
arguments are suitable for all solutions $(<?) if singularities are situated at the real axis 
of complex z-plane in particular for the iterative solutions of Refs. 0, [3|] . 

Let us see what the singularities of the integrand of (|H|) in the complex t-plane are. 
First of all the integrand has simple poles at t = Inx ± m(l + 2n), n = 0,1, 2,.... All the 
residues of the function exp(t)/(x + exp(t)) at these points are equal to unity. Apart from 
this poles the integrand of flTB] ) has at t = poles up to fourth order and logarithmic 
type branch points which coincide with poles from second order to fourth order. Let us 
cut the complex t-plane in a standard way, t = =bi7r — A, with A being the real parameter 
varying from to oo. Append the integration by the arch of the "infinite" radius without 
affecting the value of the integral. Close the integration contour C\ in the upper half-plane 
of the complex variable t excluding the singularities at t = m. In this case an additional 
contribution emerges due to the integration along the sides of the cut and around the 
singularities at t = in. The corresponding contour we denote as C 2 . 

Let us turn to the integration along the contour G\. For the integrand of Eq. ([13]) 
which we denote as F(t) the residues at t = Inx + m(l + 2n), n = 0, 1, 2, ... are as follows 

ResF(t) | 4= i na;+ i 7r (i + 2n)= $(lnx + 2mn) - $(lnx + 2m(n + 1)). (16) 



By using the residue theorem one readily obtains the contribution A(x) to the integral (|13"D 
from the integration along contour C\. It reads 

1 r 00 

A(x) = / F(t) dt=T\ Res F (t = Inx + z'tt(1 + 2n)) = $(lnx). (17) 

We can see that this contribution is exactly equal to the initial a(x). Therefore we call it 
a perturbative part of a an (x), a pt (x) = A(x). The remaining contribution of the integral 
along the contour C 2 can naturally be called a nonperturbative part of a an (x) according 
to Eq. (||). Let us turn to the calculation of d^(x). We have 

c 2 

We consider x as real variable, x > 1. One can omit the term of the integrand in Eq. ( ]I8D 
which has singularities at t — —m. Let us change the variable t = z + m and introduce 
the function 

/(*) = 1 X -^-y (19) 

1 — x exp(— z) 
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Then we can rewrite Eq. ( |T8| ) in the form 



2m 



o 



Hz) c 



r ln 2 ^ 



z 



2" 



-6 3 



ln 3 U) 



+ 3C- 



5\ In 2 



27 2 4 



+ (3C 2 - 5C + 3k + 1 



In* C 3 - \C 2 + (3k + 1)C + re' 



(20) 

The cut in the complex z-plane goes now from zero to — oo. Starting from z = — oo — iO 
the contour C goes along the lower side of the cut then goes around the origin, and then 
it goes further along the upper side of the cut to z = — oo + iO. The contour C can be 
chosen in such a way that it does not envelop a " superfluous" singularities corresponding 
to the perturbative contributions. Using the technique of Ref. @ one can integrate the 
terms of Eq. ( |20| ) and obtain 



a an ( X ) 



1 



(1 + C)x 



x — 1 



[X 



+ x [ da ln(— In a) 
1 r J 



x + a 

(x — ay 



--b 2 



7T 



X 



l-y + re + (l + C) 2 



x 2 + Axa + a 2 



x(x + 1) 



[X 



- 1 f , 

— + xj da 



2(1 + C) ln(-lna) +ln 2 (-lncr) 



(x - 



7T 



2 + -re + R + 3(1 + C) 1 - — + k + (1 + C) 



X ■ 



x(x 2 + Ax + 1 



[x 



'-f^ + xjda [s(l-^ + n + (l + C) 2 ) 



x ln(-lno-) + 3(l + C) ln 2 (- In a) + ln 3 (- In a) 



a; 3 + llx 2 a + llxa + a" 



(x - a 



:2D 



This formula gives the nonperturbative contributions in the explicit form. Expanding 
Eq. ([H]) in inverse powers of x we have 



n=l 



(22) 



Making the change of variable a = exp(— t) and integrating |L8], [19] over £ we finally find 



-1 + 6ra [1 + C - 7 - ln(n)] - ^b 2 n 2 



7T Z / 

1- — + K+ I + C- 7- ln(n 
o 



+ -b 3 n 3 
6 



2 + -re + R - 2C 3 + (l + C - 7 - ln(n)) + 3(l + C - 7 - ln(n)) ( 1 - y + k 



(23) 

Here 7 is the Euler constant, 7 ~ 0.5772. We can see from Eq. (^) that power series (|22f) is 
uniformly convergent at 2 > 1 and its convergence radius is equal to unity. For numerical 
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Table 1: The dependence of c\ on nf for 1-loop, 2- loop, 3-loop, and 4-loop cases. 



n f 


1— loop 

Cl 


A 

^2-loop 


^3— loop 


A 

4- loop 


2— loop 


3— loop 


4— loop 
Cl 





-1.0 


0.35640 


-0.01568 


-0.03900 


-0.64360 


-0.65929 


-0.69828 


3 


-1.0 


0.33405 


0.01608 


-0.07825 


-0.66595 


-0.64987 


-0.72812 


4 


-1.0 


0.31252 


0.04949 


-0.11006 


-0.68748 


-0.63799 


-0.74805 


5 


-1.0 


0.27813 


0.11653 


-0.16002 


-0.72187 


-0.60535 


-0.76537 


6 


-1.0 


0.22433 


0.25378 


-0.22731 


-0.77567 


-0.52189 


-0.74920 



evaluation of the coefficients c n assume that C — 0. Then the coefficients c n depend on 
n, rif and number of loops taken into account. The 1-loop order contribution to c n equals 
to —1 for all n and nf. Up to 4-loop approximation the coefficients c n for all n, nf are 
negative. With the exception of the 3-loop case at rif = 6, the 2 — 4-loop coefficients c n 
for rif = 0,3,4,5,6 monotonously increases in the absolute value with increasing n. In 
the ultraviolet region (x ^> 1) the nonperturbative contributions are determined by the 
first term of the series (^). In Table p] we give the values of c\ and loop corrections for 
nf = 0, 3, 4, 5, 6. One can see that for all nf up to four loops c\ is of the order of unity. 
An account for the high loop corrections results in some compensation of the 1-loop 
leading at large x term of the form 1/x. The relative error of the approximation of the 
analytic running coupling with only one first term of the series (|22|) for the nonperturbative 
contributions taken into account, 

a an {Q 2 ) * c^(Q 2 ) - ^ jl - b (1 - 7 ) + h 2 (l - ^ + k + (1 - 7 ) 2 ) 



2 + \k + k - 2C 3 + (1 - 7) 3 + 3 (1 - 7) ( 1 - + « 




(24) 



5 v " v " \ 6 

is shown in Fig. [I] for 1 — 4-loop cases. We have chosen nf = 4 within the region of Q 
considered. 

As a result of the expansion coefficients increase not too fast, the representation of the 
analytic running coupling of QCD in the form (fj), (|3]) with c n as ( |23|) provide one with 
the effective method for the calculation of a an at Q > A. At that there is no need for the 
summation of large number of terms of the series. Approximation of the nonperturbative 
"tail" by the leading term already at Q ~ 5 A gives one percent accuracy for a an . As 
regards the parameters A^ nf=5 \ A (n/=4 - 1 of the analytic approach they turn out to be close 
to the perturbative values as a consequence of the rapid decrease of the nonperturbative 
contributions in the ultraviolet region. 

I am indebted to B.A. Arbuzov, V.A. Petrov, V.E. Rochev for useful discussions. This 
work was supported in part by the Russian Foundation for Basic Research under grant 
No. 99-01-00091. 
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Figure 1: Relative error of the approximation of a an with a™^' approximated by one first 
term of the series, as function of x 1 ^ 2 = Q/A for 1 — 4- loop cases at rif — 4. Dash-dotted 
line, dotted line, dashed line and solid line correspond to 1-loop, 2-loop, 3-loop and 4-loop 
cases respectively. 
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